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ABSTRACT 



We present the first spectroscopic images of solar radio transients from 
the prototype for the Murchison Widefield Array (MWA), observed on 2010 
March 27. Our observations span the instantaneous frequency band 170.9- 
201.6 MHz. Though our observing period is characterized as a period of 'low' 
to 'medium' activity, one broadband emission feature and numerous short-lived, 
narrowband, non-thermal emission features are evident. Our data represent a 
significant advance in low radio frequency solar imaging, enabling us to follow 
the spatial, spectral, and temporal evolution of events simultaneously and in 
unprecedented detail. The rich variety of features seen here reaffirms the coro- 
nal diagnostic capability of low radio frequency emission and provides an early 
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glimpse of the nature of radio observations that will become available as the next 
generation of low frequency radio interferometers come on-line over the next few 
years. 

Subject headings: Sun: corona — Sun: radio radiation — radiation mechanisms: 
non-thermal — instrumentation: interferometers 



1. Introduction 

Low radio frequency < 300 MHz) emission provides powerful diagnostics of the so- 
lar corona. However, high fidelity solar imaging at low frequencies is challenging. Coronal 
emission features are complex and dynamic, evolving rapidly in space, time, and frequency. 
Consequently, the limited instantaneous spatial and frequency coverage provided by cur- 
rent radio interferometers have been inadequate to simultaneously resolve transient solar 
phenomena spatially, temporally, and spectrally. 

This situation should change dramatically in the next few years as a new generation 
of low-frequency radio arrays becomes available, leveraging recent advances in digital signal 
processing; hardware a nd computational capacity. The Murchison Widefield Array (MWA) 



(iLonsdale et al.l 120091 ) will be one such array, with most of its 512 elements spread over 
1.5 km and a few outliers out to 3 km. The resulting dense instantaneous monochromatic uv 
coverage will provide a radio imaging capability with unprecedented fidelity and fiexibility. 
The MWA is currently under construction at the Murchison Radioastronomy Observatory, in 
the remote and radio-quiet Western Australian outback, and a prototype system comprising 
32 interferometer elements (tiles) is operational on site. This system, hereafter referred to 
as the "32T" , serves as an engineering testbed and provides early science opportunities, in 
advance of the full MWA. We present here solar imaging observations obtained with the 32T. 
The quiescent solar emission in the MWA frequency range of 80^300 MHz is dominated by 
coronal emission from heights of ^HO above the photosphere. These are among the 
first high-fidelity, high dynamic-range, spectroscopic images of the Sun with a good temporal 
and spectral resolution at meter wavelengths. 



2. Data 

The 32T tiles are arranged along the arms of a randomized Reuleaux triangle (see 
Cohanim et al. 2004), providing a fairly uniformly sampled uv plane with baseline lengths 
up to 350 m (Fig. [1]). For solar observations, the correlated fiux is dominated by the Sun, 
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reducing the effective field-of-view (FoV) to ^1°, for which the uv samphng provided by the 
32T exceeds the Nyquist crite rioijll. This, together with a high signal-to-noise ratio, permits 



robust, high fidehty imaging ( IBraceweU fc Robertslll954l ). 



In its usual interferometric mode of operation, the 32T provides auto- and cross-correlations 
(visibilities) for all 64 input signals (32 tiles x 2 linear polarizations). The entire RF band 
is directly sampled and filtered into 1.28 MHz wide coarse channels. Twenty four of these 
coarse channels are further processed by the correlator, which can currently provide a time 
resolution of 50 ms with a 50% duty cycle over a 30.72 MHz band and a spe ctral resolution 



of 40 k Hz. The 32T design closely follows the MWA architecture described in lLonsdale et al 



(120091 ). 



The data presented here were obtained on 2010 March 27 from 04:24:53 to 04:34:03 UT 
and span the range 170.8-201.6 MHz. They were smoothed to a time resolution of 1 sec- 
ond before further processing. Data editing, calibration, and imaging were performed using 
the Astronomical Image Processing System. Self-calibration, to solve for frequency inde- 
pendent complex gains, was performed individually for each of the coarse channels using a 
10-second interval where the Sun was in a relatively quiescent state. These gain solutions 
were then applied to the full time interval. No absolute fiux calibration was possible, owing 
to the lack of observations of a suitable calibrator. The amplitude part of the bandpass 
was calibrated using the total power spectra from the same quiescent time interval as used 
for self-calibration. As will become evident later, the choice of a spectral index only im- 
pacts the underlying spectral slope. We assumed a spectral index a=2.6 for the quiet Sun 
{Siy oc jy^ where 5^^^ is the fiux density and u is the observing frequency); this choice is based 
on 32T observations during the rece nt deep solar mi nimum (2008 November), which yielded 



a = 2.6±0.4, a value consistent with lErickson et al.l (119771 ). All the data presented here cor 



respond to the east-west (XX) polarization. The wide FoV MWA tiles are expected to have 
significant (but stable) cross-polarization leakage, and an absolute polarization calibration 
was not attempted. 

Imaging and deconvolution were performed using the standard CLEAN algorithm with 
robust (7^=0) weighting, resulting in a synthesized beam at the band center of 961^' x 796^' 
(FWHM). For the images presented here, the data were averaged over a single coarse channel 
during the gridding process and restored with a circular beam with FWHM 800'^ The edges 
of the coarse channels could not be calibrated satisfactorily, hence the first six and the last 
four spectral channels of each coarse channel were fiagged during imaging. 



^Su ^ 1/2^07 where Su is length scale in the uv plane beyond which visibilities are no longer correlated, 
and ^0 is the size of the FoV in radians 
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Fig. 1. — The instantaneous uv coverage for a single coarse channel centered at 186.2 MHz. 
The circle indicates the size of a uv-cell^ corresponding to the Nyquist sampling criterion for 
a 1 degree FoV. 
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Fig. 2. — Observed visibilities as a function of frequency (XX polarization). The left panels show amplitude (arbitrary 
units), and the right panels phase (in degrees). The top and bottom rows correspond to baselines with projected 
lengths of ^40A and ^196A, respectively, at the band center (186.26 MHz). The frequency (x-axis) ranges from 170.8 
to 201.6 MHz; the time (y-axis) spans 550 s starting at 04:24:53 UT on 2010 March 27. The horizontal bar below each 
panel shows the color coded scale and the dark vertical streaks arise from flagging of the coarse channel edges (Sec. [2j). 
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3. Analysis and Results 

Amplitudes and phases observed on a representative short (^64 m) and long (^315 m) 
projected baseline are shown in Fig. [21 Baselines of different lengths and orientations measure 
independent Fourier components of the source structure. The examples shown here illustrate 
several key features of the data. Only one polarization is shown as both polarizations show 
very similar behavior. The most prominent feature in Fig. [21 visually obvious on all baselines 
in both amplitude and phase, spans the entire observing band around 04:30:10 UT. We 
hereafter refer to this as the "broadband" feature. Its properties are similar to a weak 
type III burst (Sec. [Ij). In addition, a large number of shorter-lived, narrowband features are 
evident in the visibility amplitudes. Many of these amplitude variations are accompanied 
by corresponding variations in phase, implying a change in the brightness temperature (T^) 
distribution in the corona. These features become more numerous and prominent visually 
with increasing baseline length and are referred to as "narrowband" features hereafter. A 
dearth of truly "quiescent" periods is also evident, and the ubiquitous modest variations 
seen in the spectral structure are solar, not instrumental, in origin. 

To further illustrate the diversity in the types of emission seen in the short span of data 
presented here, as well as its spectrally complex and highly time-variable behavior. Fig. [3] 
shows several series of autocorrelation (total power) spectra. The bottom panel shows one of 
the few quiescent intervals; the baseline spectral slope seen here is due to the spectral index 
of the quiet Sun, dominated by thermal emission. An amplitude variation of ^5% over a 10 
s interval is seen even in "featureless" parts of the data. The middle panel shows a series 
of spectra illustrating one of the numerous "narrowband" features seen in Fig. [2 These 
features show remarkably rapid evolution in spectral shape and intensity. They typically 
span ^5-10 MHz in bandwidth, outside of which the spectral flux density returns to that 
of the quiet Sun. Their peak flux density approaches ^1.5 times the quiescent solar flux 
density at that frequency. Finally, the top panel of Fig. [3] shows the dynamic behavior of the 
"broadband" feature around 04:30:10 UT. Its bandwidth exceeds the ^30 MHz observing 
bandwidth, and at its peak, its flux density is ^2.5 times the quiescent solar flux density. 
The abrupt changes from one second to the next imply that the physical changes leading to 
the production of this emission are temporally undersampled. 

The data presented here permit a high fldelity image of the Sun for every individual time 
and frequency slice (i.e. every pixel on the dynamic spectra shown in Fig. [2]). Fig. [H shows 
a collage of images that illustrates the rapid evolution of the appearance of the solar corona 
during our observations. Although the spatial resolution of the 32T is limited, the fldelity 
and dynamic range of these images is unprecedented at these frequencies; the dynamic ran^e 



of ^2500 exceeds that of earlier images by about an order of magnitude (cf. iMercier et al. 
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Fig. 3. — Sample autocorrelation spectra from tile 6, X polarization. Top panel: a 14 s 
period bracketing the most prominent "broadband" feature seen in Fig. [2l Middle panel: 
variations seen across 10 s for one of the numerous shorter-lived "narrowband" intensity 
enhancements. Bottom panel: a 10 second interval exhibiting little intensity variation. The 
inset in the bottom panel zooms in on a part of the band. As with Fig. O the data gaps 
reflect flagged data. 
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Fig. 4. — Top and middle rows: sets of images, 1 s apart, from the "broadband" and 
one of the numerous "narrowband" features, respectively, highhghting the rapid temporal 
variation in the emission morphology. Animations including these panels are available in 
the on-line journal. Bottom row: A 32T radio image of the quiescent Sun (left) and a 304A 
image from SOHO/EIT (right), taken a few hours earlier (01:19 UT). The bright region in 
the northeastern quadrant is the region 11057. The red circle on the 32T image represents 
the size of the optical solar disc and has been registered by aligning the centroid of radio 
emission at the location of the active region 11057 with its counterpart on the 304A image, 
ignoring the time offset. All radio images are at 193.3 MHz (XX polarization) and have the 
same intensity scale (0-2500, arbitrary units). All images have celestial north on top. 
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20061 : iMercier fc Chambd l2009l ) . The higher fidehty and dynamic range enables us not only 
to image the quiescent sol ar emission in the presence of niuch brighter feature s , unlike most 



earlier observations (e.g., iKail Il970l : iKundu et al.l 1 19861 : iKrucker et al.l 1 19951 : IVilmer et al. 



2OO2I : iRamesh et al.ll2010l ). but also to track low level variations in the coronal emission. 
The bottom panels in Fig. H] illustrate the morphological relationship between our 32T radio 
images and the extreme ultraviolet (EUV) images taken on the same day. A persistent feature 
coinciding with the location of the NOAA Space Weather Prediction Center (SWPC) region 
11057 (hereafter region 11057), is seen in the solar radio images, even at quiescent times. 

Rapid and significant changes in the radio morphology of the Sun on timescales of 
seconds are clearly visible in the images presented in Fig. d] and are further highlighted in 
the two "movies" available in the on-line journal. A significant fraction of this variability 
maps back to the location of the region 11057. Subtle changes in Tb as a function of time 
are also discernible at the locations of the active regions on the Eastern limb and close to 
the center of the disc. 

Fig. [5] highlights the spectroscopic capability of the 32T by showing the spectra from 
representative pixels across the entire emitting region. Various interesting features can be 
seen here. Even at quiescent times, there is a significant difference in the spectral slope 
of regions close to and away from region 11057. For both the "narrowband" and "broad- 
band" features, the changes in spectral structure are confined to the vicinity of region 11057, 
and to a smaller extent to the locations of other active regions visible in the EUV image 
(Fig.Hj). The "narrowband" features are consistent with spectral features superposed on the 
underlying quiescent spectrum. In contrast, the underlying quiescent spectrum is no longer 
visible in spectra of the "broadband" feature, and the changes in spectral characteristics are 
significantly larger, more dynamic, and complex. 



4. Discussion and Conclusions 



At low radio frequencies, a featureless continuum due to thermal bremsstrahlung from a 
corona with optical depth approaching unity cons titutes the bulk of the "qui et" Sun emission 
and is expected to vary only slowly in time (e.g., ISheridan fc McLeanlll985 l). In this regime, 
emission that exhibits rapid spectral variations is dominated by plasma emission, a resonant 
emission mechanism in which electrostatic Langmuir waves are converte d to electromagnetic 
radiation at the local plasma frequency, /p, and its harmonic, 2/p (e.g.. iRobinson fc Cairns 



2000). This is the accepted er aission mechanism fo r type II (e.g., iNelson fc Melrosdll985l ) 



and type III solar bursts (e.g.. ISuzuki fc Dulklll985l ) 
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Fig. 5. — Spatially localized spectra of emission across the solar disk at several of the same 
times shown in Fig. [H Celestial north is on top, and the red circle marks the size of the 
optical solar disc. Spectra are shown for every third pixel, where the pixel size is 100''xlOO'^ 
The spectra span the entire observing bandwidth, binned into 24 frequency bands spaced 
1.28 MHz apart and averaged over 0.8 MHz. The y axis ranges are (in arbitrary units): top 
row: 50-4996 and 50-4620, respectively; middle row: 100-6000 for both panels; bottom row: 
100-1800 for both panels. Animations showing additional frames are available in the on-line 
journal. The frames corresponding to the bottom two rows have fixed ^-axis scales; those 
for the top row are autoscaled. 
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The ability of the 32T to image and track the temporal and spectral evolution of coronal 
emission with high dynamic range and fidelity over a broad spectral band represents a major 
advance in the quality, and consequently the utility, of radio measurements for understanding 
the dynamic solar corona. For example, we see numerous instances of transient, non-thermal 
emission in our short observations, many of which outshine the Sun by a factor of few and 
dramatically alter the Tb distribution. Though there have been earlier imaging observations 
of solar radio tranisents, these investigations have been limited to a handful of well separated, 
non-contiguous frequencies due to the limited spectral sampling provided by previous radio 
interferometers. 

The coherent emission mechanisms inherent at low radio frequencies have long been 
known to be a sensitive probe of electron acceleration in the the corona. Given our much 
improved observing capability, it is not a surprise that in spite of significant transient be- 
havior seen in our data, the level of solar activity on the day of our observations has been 
characterized as /owl to medmTril. The SWPC "events list" reports nine type III radio bursts 
on this day, none of them during the interval presented here. The SWPC also reports seven 
GOES B-class and four GOES C-class X-fiares on this day (ten and five X-ray fiares, respec- 
tively in a 24 hour period around our observations). All of the fiares for which positional 
information is available (C-class) were associated with the region 11057. The start time for 
one of the B-class fiares coincides with the stop time of our observations. 

Based on its appearance in the frequency-time plane, the "broadband" feature present 
in our data (Fig. [2j) r esembles a weak typ e III burst and shares similarities with the "mi- 
crobursts" reported bv lKundu et al.l (119861 ). A corresponding feature was seen at low signal- 
to-noise levels between ^80 and 150 MHz by the Learmonth and Culgoora radiospectro- 
graphs, suggesting that this feature peaked in intensity within that spectral window. Though 
it appears to be a feature with a fast spectral drift rate, the signal-to-noise ratio and tem- 
poral sampling of the Learmonth and Culgoora data are not sufiicient for a reliable estimate 
of the drift rate. Nancay Radioheliograph data spanning the wide 150^455 MHz frequency 
band show two persistent but time variable features separated by ^ 200^' associated with 
region 11057 at nearby times. This suggests the presence of an underlying weak type IV 
emission. The non-thermal spectra of the region 11057 even at "quiescent" times, and the 
persistent feature seen at this location in our data are consistent with this. 

Lacking an absolute fiux calibration, we compute an approximate Tg of the wide- 
band burst at the center of our observing band by assigning the "quiet" Sun a Tb of 
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?^7.5xl0^±15% K, based on measurements of iLantos fc AvignonI (119751 ) at 169 MHz, and 
scaling it to the peak of the burst emission (04:30:12.5 UT). A measured receiver tempera- 
ture of ^50 K and an estimated Galactic background contribution of ^150 K were used. We 
estimate an average Tg ^3.1x10^ K, signi ficantly lower than th e average value of 6.3x10'' K 
reported for type III bursts at 169 MHz (jSuzuki &: Dulklll985l ). The low value of Tg mea- 
sured highlights the ability of the 32T to make sensitive measurements. We note that since 
the angular size of the emitting region is expected to be much smaller than the 32T beam, 
the measured Tb is expected to be considerably reduced due to beam-dilution. Assuming a 
plasma emission process and the Newkirk (lOxNewkirk) model (jNewkirklll96ll ) for coronal 
electron density leads to a coronal height of 1.7 Rq (2.8 Rq) above the photosphere for the 
source region at our central observing frequency. A bandwidth of 30 MHz corresponds to a 
range of ^0.09 (^0.25 R©) around this coronal height. 

A particularly intriguing aspect of our data is the abundance of short-lived, "narrow- 
band" emission features carpeting the frequency-time plane, which to our knowledge, have 
not been reported previously. Though smaller in spectral and temporal extent, these features 
have a similar Tg, implying a coherent emission process, like plasma emission. The existence 
of these emission features provides evidence for impulsive non-thermal energy deposition 
into the corona at levels lower than previously known. Evaluating their contribution to the 
coronal heating budget will require building up a better understanding of the frequency of 
occurrence and energetics of such features. Existing radio observations have been shown to 



be a more sensitive probe for such phenomena than X-ray observations (e.g. iRamesh et al. 



20101 ). and these observations extend the radio sensitivity advantage by another order of 
magnitude. 

As seen in Fig. HJ much of the "narrowband" activity seen in Fig. [2] maps to the lo- 
cation of the region 11057. Even the low level fiuctuations in the spectrum when the Sun 
is in a relatively quiescent state (e.g.. Fig. [3l bottom panel; Fig. [5j), manifest themselves 
in the image plane as flickering intensity at the location of this active region. The occur- 
rence of many weak X-ray flares from region 11057 around our observing period suggests 
a persistent low-level activity during our observations. This raises the possibility that the 
"narrowband" features are the radio signatures of the same magnetic fleld rearrangement 
activity: numerous small reconnection events too faint to picked up by X-ray monitoring 
platforms like GOES, but with weak type Ill-like radio emission strong enough to be seen by 
the 32T. This hypothesis is strengthened by the observation that the "narrowband" features 
were absent on 2010 March 29 (not shown here), a day when no X-ray activity was reported. 
There have been attempts to investigate the possibility of usinR type III solar bursts as pre- 
dictors for large solar flares in the past (j Jackson fc Sheridanlll979l : lKandll98ll ). The ability 
to characterize type Ill-like activity in much greater detail and with much higher sensitivity. 
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demonstrated here, enhances the prospects for investigating such a connection. 

Compared to the 32T, the MWA will provide an order of magnitude improvement in 
angular resolution, 16 times as much collecting area and 16^ times as many baselines. With its 
unprecedented high fidelity and dynamic range polarimetric spectroscopic imaging capability 
for every spectral and temporal slice, and a radio quiet site providing seamless access to the 
RF spectrum, the MWA will be an excellent match to the challenge of solar imaging. Every 
significant advance in our capability to observe the Sun has revealed that the Sun and 
corona are more dynamic and feature-rich than had been thought previously. HINODE and 
the Solar Dynamics Observatory (SDO) represent the latest in that progression. The 32T 
array is already demonstrating that we are beginning to see features at a level of detail not 
seen before, and we can justifiably expect the MWA to similarly add to our understanding 
of the solar corona and the heliosphere. 
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AST-0457585 and PHY-0835713), the Austrahan Research Council (grants LE0775621 and 
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